Limits on diffuse fluxes of high energy extraterrestrial neutrinos with the 

AMANDA-BIO detector 
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Data from the AMANDA-BIO detector taken during the austral winter of 1997 have been searched 
for a diffuse flux of high energy extraterrestrial muon-neutrinos, as predicted from, e.g., the sum 
of all active galaxies in the universe. This search yielded no excess events above those expected 
from the background atmospheric neutrinos, leading to upper limits on the extraterrestrial neutrino 
flux. For an assumed spectrum, a 90% classical confidence level upper limit has been placed at 
a level E^"1>(E) — 8.4 x 10~^ cm"'^ s~^ sr~^GeV, (for a predominant neutrino energy range 6-1000 
TeV) which is the most restrictive bound placed by any neutrino detector. When specific predicted 
spectral forms are considered, it is found that some are excluded. 

PACS numbers: 95.55.Vj,95.85.Ry,96.40.Tv,98.54.-h 



High energy extraterrestrial neutrinos are believed 
to be produced in energetic accelerated environments 
through proton-proton or proton-photon interactions via 
pion production and decay. Such an accelerator might 
be the core of an active galaxy, powered by a supermas- 
sive black hole. In theirpioneering work, Stecker, Done, 
Salamon and Sommers calculated the expected diffuse 
flux of neutrinos from the sum of all active galaxies and 



found that such a flux could be observable deep under- 
ground in a large neutrino detector. Further predictions 
have followed (for a summary see for example the review 
of Learned and Mannheim and with the construc- 
tion and operation of the first high energy neutrino de- 
tectors, the sensitivity has been reached to enable such 
predictions to be tested. Searches have been made and 
limits have been reported by the DUMAND Q, Frejus 
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0, Baikal(i/e) 0,11], MACRO and AMANDA(r/e) 
neutrino detectors. In this letter, we describe the search 
for high energy extraterrestrial neutrino-induced muons, 
using data collected during the austral winter of 1997 
with the AMANDA-BIO detector j9|,Uj|, located in the 
antarctic ice cap at the South Pole station. 

The AMANDA (Antarctic Muon And Neutrino Detec- 
tor Array) telescope detects high-energy muon-neutrinos 
by observing Cherenkov light from muons resulting from 
neutrino interactions in the ice surrounding, or the rock, 
below the detector. While extraterrestrial neutrinos will 
produce high-energy muons from all arrival directions, 
those coming from above the detector will be very diffi- 
cult to separate from the overwhelming flux of downward- 
going cosmic-ray induced atmospheric muons. The ma- 
jority of these muons are rejected by only accepting up- 
ward going neutrino-induced muons; the earth filters out 
muons produced in the atmosphere on the other side 
of the planet. There is a small remaining flux of mis- 
reconstructed events which is removed by quality cuts 
that leave only well reconstructed events. After the at- 
mospheric muons are removed, there will remain a flux 
of upward going muons from cosmic ray induced atmo- 
spheric neutrinos that have penetrated the earth and in- 
teracted near the detector. This relatively well under- 
stood neutrino flux is a background to the search but 
has been used to verify the performance of the detector 
0, 0| . The separation of the extraterrestrial neutrino- 
induced muons from the atmospheric neutrino-induced 
muons is based on the expected energy spectrum of the 
detected muons. Typically, a model of an extraterrestrial 
source of neutrinos has a harder spectrum (e.g. E^^) 
111 than that of the atmospheric neutrinos (^ E^"^-^) 
12. Ill^ . After accounting for neutrino interaction and 
muon propagation, this energy difference carries over to 
produce a harder muon energy spectrum for the extrater- 
restrial neutrino-induced muons near the detector. The 
energy of the muon is not measured directly, but more 
energetic muons tend to produce more Cherenkov light 
and thus more hit optical modules in the detector, and 
this observable, the channel multiplicity (Nch), is used as 
the primary separator of higher energy extraterrestrial 
neutrino-induced muon events from the background of 
lower energy atmospheric neutrino-induced muons. 

In this analysis, the event selection cuts were designed 
to retain high energy track- like events Q|. The detec- 
tor simulation has changed from that used in the at- 
mospheric neutrino analysis |lfll |. A new muon propa- 
gation code was used, which accounts for all rel- 
evant stochastic light emission from the muons. The 
depth-dependent optical properties of the fiducial ice 
were determined using atmospheric muons as a calibra- 
tion source. 

Before the energy sensitive channel multiplicity cut was 
finally applied, 69 events remained in the data sample, 
whereas a full simulation of the detector response to the 
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FIG. 1: Channel multiplicity distribution after final cuts, 
showing the expected excess of events from an E"'^ spectrum 
at the higher multiplicities. 



atmospheric neutrino (Lipari |l2j') flux (neglecting neu- 
trino oscillations, which would reduce the prediction by 
only a few percent) predicts 85 events for the 130 days 
of live time. The absolute difference in the numbers of 
events is consistent with Poisson fluctuations, or with the 
±25% [l^ uncertainty in the atmospheric neutrino flux, 
or with uncertainties in the simulation efficiencies. The 
distribution of the data and atmospheric simulation are 
shown in figure ^ The error bars on the data are 90% 
unified confidence intervals [l^ for the fixed but unknown 
value of the mean rate = ^{,-|-/is for each bin. Only one 
bin (Nch = 25 — 30) has a background prediction incon- 
sistent with the confidence interval. More specificially, a 
generalised likelihood ratio test of the shape of the at- 
mospheric neutrino hypothesis as the parent distribution 
of the data yields a chance probability of 20%, which 
is too large to reject the shape of the atmospheric neu- 
trino hypothesis. We choose to treat the rate of observed 
atmospheric neutrinos as a constraint on the overall de- 
tector efficiency and then carry through an efficiency un- 
certainty from the atmospheric neutrino flux prediction 
and Poisson error on the observed rate. Therefore, to 
calibrate the overall detector sensitivity, we take the 69 
events as the best-fit estimate of the number of atmo- 
spheric neutrinos and rescale all efficiencies by a factor 
69/85. This is conservative, since if the first bin dis- 
crepancy was due e.g. to a simulation effect, then no 
renormalisation would be needed, and the limits would 
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improve slightly. We combine the Poisson error on the 
observed rate of atmospheric neutrinos with the theo- 
retical flux uncertainty (taken as a uniform probability 
distribution centred about the best-fit flux $ and extend- 
ing to ±0.25$) to compute the correlations between the 
background and efficiency for later use in the probability 
distribution function used in the confidence interval con- 
struction. To incorporate these systematic uncertainties 
in the efficiencies into the limit calculati ons, we follow 
the prescription of Cousins and Highland as imple- 
mented by Conrad et al. [l^ with the unified Feldman- 
Cousins ordering and improved by a more appropriate 
choice of the likelihood ratio test We also report 

all limits and sensitivities with and without the assumed 
uncertainty. 

In addition to the data and atmospheric neutrino pre- 
diction, figure n also shows the prediction for an 
signal flux at a level $(E) 10^^ cni"^ s"^ sr~^ GeV, 
a flux that would have been readily detected. Setting a 
limit on a flux <i>(E) involves determining an experimental 
signal event upper limit /i(nobs,nb), which is a function 
of the number of observed events, Uobs) and expected 
background, Ub, after the cuts are applied. A simula- 
tion chain accounting for neutrino absorption, interac- 
tion and neutral current regeneration, muon propagation 
and detector response gives the number of signal events. 
Us, expected from the source flux <I>(E). The limit on the 
source flux will then be $iimit (E) = (E) x (uobs , nb ) /ug . 
The choice of final cut is optimised before examining the 
data by minimising the average "model rejection factor" 
(MRF) /l(nobs, nb)/ns where the as yet unknown 

experimental event limit /zfunKs, Uh) is replaced by the 
average upper limit /2(nb) [ig- Over an ensemble of hy- 
pothetical repetitions of the experiment, this choice of 
cut will lead to the best average limit $iiniit(E). 

The integrated channel multiplicity distribution is 
shown in figure |21 Also shown is the 90% confidence 
level Feldman- Cousins average upper limit which is a 
function of the expected background. The optimal cut 
is the one where the model rejection factor /i(nb)/ns is 
minimised. Figure |21 also shows the average flux upper 
limit (E^$ X MRF) as a function of the choice of mul- 
tiplicity cut. The minimum flux limit occurs at a cut 
of Nch > 54, where we expect nb — 3.06 and an av- 
erage signal event upper limit of 4.43 ignoring the un- 
certainties in the efficiency and background, and 4.93 
when the uncertainties are included. The 10~^E~^ sig- 
nal flux would produce 56.7 events. This leads to cor- 
responding expected average limits on the source flux of 
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E2$9o%(E) = 7.8 X 
uncertainties), and 1 
ing uncertainties). 

We note that the expected overall flux limit is rela- 
tively insensitive to the choice of cut, with a broad mini- 
mum seen in flgure[21in the range of multiplicities 50-70. 
We now apply this optimal multiplicity cut to the data. 
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FIG. 2: Integrated distributions of event numbers as a func- 
tion of the multiplicity cut (top plot). The minimum in the 
flux average upper limit (bottom) is found by minimising the 
ratio of the average event upper limit to the expected E"^ 
signal. 



and flnd that 3 events remain. Ignoring the systematic 
uncertainties gives an event limit of 4.36 and a flux up- 



7.7 X lO-'^cm-^s" 



sr-^GeV. 



per limit of E2$9o%(E) 
Including the systematic uncertainties leads to an event 
limit of 4.75 and our flnal flux limit on an E"^ spectrum 
of E2$9o%(E) = 8.4 X 10-^cm-2s-isr-iGeV. 

Figure |31 shows the neutrino energy spectrum of the 
simulated events before and after the multiplicity cut of 
54 channels, for both atmospheric neutrinos and neutri- 
nos from an E^^ spectrum. The multiplicity cut corre- 
sponds to a sensitive energy range of 6-1000 TeV, which 
contains 90% of the expected E"^ signal. The peak re- 
sponse energy is just below 100 TeV. 

Just as a limit was placed on an assumed E"^ spec- 
trum, limits can be placed on any neutrino flux pre- 
diction, and we consider a sample of predictions that 
are near the limit-setting capability of the detector 
(MRF '^1). For each case, we optimise the flnal Nch cut 
to minimise the expected average flux upper limit, then 
compare the expected number of extraterrestrial neutrino 
events after the cuts to the observed event limit; those 
predictions that produce expected event numbers greater 
than the observed event limit are excluded at the stated 
classical confidence level. The results of these calcula- 
tions are shown in tableland in figure 01 For each flux 
we again report two sensitivities and limits - one assum- 
ing no systematic uncertainties and the second including 
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FIG. 3: Energy spectrum of the incident atmospheric (dashed 
lines) and (solid lines) neutrinos for events that pass the 
initial cuts (thin lines), and have channel multiplicity greater 
than the optimum cut of 54 channels (thick lines) . 



systematic uncertainties. We find that the predictions of 
Szabo and Protheroe (SPH92L ,21|], P96p7pp ^) are 
excluded. The quasar core (SSQC) prediction of Stecker 
and Salamon is just excluded (MRF = 0.98), but 
the blazar jet (SSBJ) prediction is not. The limit of the 
original Stecker, Done, Salamon and Sommers flux 1] 
(SDSS) is a factor of 2 above the prediction and there- 
fore the prediction is not excluded. 

Maximal mixing from neutrino oscillations [23I be- 
tween and Vt during propagation to the earth would 
lead to limits a factor of two greater than quoted here 
due to the loss of half the expected flux into the i^r 
state. However, some of these Vt would regenerate in 
the earth (i^r ^ t ^ v^) lessening the effect H^H^. 

We also place a limit on a model of prompt charm in- 
duced neutrinos [l^l (ZHV92) in the earth's atmosphere 
and find that the detector sensitivity is about a factor 
of 4 away from excluding the prediction. More recent 
predictions are even further below the sensitivity of the 
detector |2^ . 

Since most events will originate from neutrinos near 
the peak of the detector sensitivity Yi^ ^ lif' GeV, the 
limits at that point for the three different spectral shapes 
(E^^, SSQC and Charm D) are similar, as seen in 
figure 01 

The limits presented in this letter, based on the first 
real-time year of operation of the AMANDA-BIO de- 



FIG. 4: Summary of experimental 90% classical confidence 
level flux limits from various detectors assuming an E~'^ spec- 
trum. From top: AMANDA-BIO [ue) Frejus 0, MACRO 
1% Baikal 6] and AMANDA-BIO (!/p)"(this work). The back- 
ground atmospheric neutrinos 12] are indicated by the hashed 
region representing the angular dependence of the flux. Also 
shown are the predicted fluxes (dashed), and AMANDA-BIO 
experimental flux limits (solid) for a difi'use AGN prediction 
(SSQC |Tll | - nearly overlapping dotted and dashed curves - 
MRF = 0.98) and for one prediction of prompt charm neu- 
trino production in the earth's atmosphere [27f . Since most 
events will originate from neutrinos near the peak of the de- 
tector sensitivity (E^/ ~ 10^ GeV), the limits at that point for 
difi'erent spectral shapes are similar. 



tector, are the strongest placed to date on extraterres- 
trial diffuse neutrino fluxes. Since that year, we estimate 
that about 10 times the exposure has been achieved in 
total with AMANDA-BIO (1997-99) and the expanded 
AMANDA-II detector (2000-present). We anticipate this 
combined data set to have a limit-setting potential more 
than three times better than the results presented here. 
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TABLE L Flux limits calculated for individual models of diffuse neutrino emission. The optimal Nch cut, expected background 
and signal for each model are shown. The average upper limit (/i(nb)) and average model rejection factor (/i(nb)/ns) are shown 
with and without the inclusion of systematic uncertainties. Finally the experimental limits (event limit = ^(no,nb)) and 
model rejection factor (/io/us) are given for both systematic uncertainty assumptions. 
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